Metastability and paramagnetism in superconducting mesoscopic disks 
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A projected order parameter is used to calculate, not only local minima of the Ginzburg-Landau 
energy functional, but also saddle points or energy barriers responsible for the metastabilities ob- 
served in superconducting mesoscopic disks (Geim et al. Nature 396, 144 (1998)). We calculate the 
local minima magnetization and find the energetic instability points between vortex configurations 
with diff'erent vorticity. We also find that, for any vorticity, the supercurrent can reverse its flow 
direction on decreasing the magnetic field before one vortex can escape. 

PACS numbers: 74, 74.60.Ec, 74.76.-w 



o 
o 

Oh' 



I 

O 

o 



> 

m 

00 

o 

OS 

-I— > 



I 

O 

o 



X 



The interest in understanding the creation and anni- 
hilation mechanisms and, in general, the stability of vor- 
tices in superfluids has been recently boosted by a series 
of technological advances in both mesoscopic supercon- 
ductors ||-||] and atomic condensates ||]. Most of the 
proposals for the creation of vortices in atomic conden- 
sates face, at the present time, severe technological dif- 
ficulties. On the contrary, mesoscopic superconductors 
in magnetic fields are already proving to be an ideal sce- 
nario where the detection and even manipulation of vor- 
tices at the individual level is becoming more and more 
feasible A good example, although not the only 

one, of single-vortex sensitivity can be found in the Hall 
magnetometry measurements performed on mesoscopic 
Al disks by Geim et al. [||]. Both field-cooled (FC) and 
constant temperature (CT) magnetization measurements 
provide evidence of the quantization of the vorticity of 
the order parameter. When the system is kept out of 
equilibrium, it can show paramagnetic response both in 
the CT and FC cases, whereas, as expected, equilibrium 
measurements always exhibit diamagnetism. 

Geim's experiment sheds light on the controversial 
paramagnetic Meissner effect (PME) confirming that the 
PME, at least in mesoscopic samples, is linked to a non- 
equilibrium magnetic response of the system. In addition 
to this, it raises a number of fundamental questions con- 
cerning the metastability of few-vortex states in meso- 
scopic systems: (i) What is the origin of the metasta- 
bility of a given vortex state?, (ii) why can mctastable 
states exhibit paramagnetism?, and (iii) what is the na- 
ture of the instability that triggers a change in the num- 
ber of vortices?. We anticipate the answers to these ques- 
tions, (i) The metastability results from the sample sur- 
face which favors a large surface superconductivity and 
opposes both vortex escape and entrance. This translates 
into a very stable vorticity or topological charge, Q, asso- 
ciated with all the local minima of the energy functional, 
(ii) On decreasing the magnetic field, a reversal in the 
direction of the total supercurrent flow associated with 



most local minima can take place before vortices escape. 
In general, and in contrast with recent work, detector 
effects need not be invoked to explain the paramag- 
netic response, (iii) The ultimate mechanism responsible 
for an instability in the vorticity or change in the number 
of vortices of the system is the disappearance of the sad- 
dle point separating a local minimum from a neighboring 
one with different Q. This is called energetic instabil- 
ity. We argue below that, either increasing or decreas- 
ing the magnetic field, there seems to be experimental 
evidence that this energetic instability is preempted by 
other mechanism, either associated with noise or some 
other type of fiuctuations (thermal or quantum). 

We start from the Ginzburg-Landau functional for the 
Gibbs free energy difference between the normal and su- 
perconducting states in an external magnetic field H: 
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where ^'(r) is the order parameter or Cooper pair wave 
fimction, /i(r) = V x A(r), and a and j3 are the conden- 
sation and interaction energy parameters, respectively. 
Numerical minimization procedures have been used in 
the past |^-|| to find global and even local minima of 
the Ginzburg-Landau functional applied to mesoscopic 
superconducting disks. However, saddle points or en- 
ergy barriers, which are essential for the analysis of the 
stability of the local minima, cannot be easily obtained 
from these methods. Before going into the details of how 
to overcome this problem, a few comments are in or- 
der. It is well known that the magnetic response to an 
external magnetic field of a superconductor varies with 
its size, geometry and orientation with respect to the 
field in a non-trivial way. Type-I superconductors expel 
the magnetic field below the critical temperature, but, 
for the largest disks in Geim's experiment, the magnetic 
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field can penetrate the interior |^,^. Moreover, multiple- 
vortex structures (7[|^, only expected in Type-II super- 
conductors, clearly reflect in the FC measurements ^] 
(see below). In summary, these disks behave like Type-II 
superconductors, which makes it possible, in a range of 
fields, to consider a uniform magnetic induction h{r) = B 
Next, we project the order parameter onto the "low- 
est Landau level" subspace [p|,pO[: 
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This subspace is spanned by normalized eigenfunctions 
of the linearized differential Ginzburg-Landau equations 
where L is the angular momentum (> 0) and ^^(r) is 
the associated nodeless function subject to the boundary 
conditions of zero current through the surface. In this ex- 
pansion Cl = ICile**^^ are complex coefficients and the 
radial unit is the magnetic length ^ — \fe*hJcB. We are 
considering the thickness of the disk to be smaller than 
the coherence length so that the system becomes effec- 
tively two-dimensional. This expansion has been shown 
to give good qualitative as well as quantitative results for 
the equilibrium properties at moderately high fields 

The central idea in our method for finding generic 
stationary solutions of the Ginzburg-Landau functional 
is to project the order parameter onto smaller sub- 
spaces spanned by a finite number N of eigenfunctions, 
{Li, L2, . . . , Lat}, where < Li < L2 < . . . < Xat. We 
will restrict our discussion to a disk radius R = 5^(0) 
which approximately corresponds to the largest disk in 
Geim's experiment |3j [^(0) is the coherence length at 
T = 0]. For such a disk size subspaces with dimen- 
sion iV > 3 do not play any role Q and the projected 
Ginzburg-Landau functional reduces to 
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where the Gibbs free energy is expressed in units of 
iJ^2^/87r {V being the volume of the disk), tL{B) is 
the energy of the quantum state L expressed in units 
of TiWc/S {ijJc — e*B /m*c)^ R is expressed in units of the 
coherence length ^ (T) and B and H are given in units of 
Hc2 [T] ■ The interactions appear in II (B) = J dr r $| , 
which can be interpreted as the interaction between 
Cooper pairs occupying the same quantum state L, and 
lLiLj{B) = / dr r $1 $1 . , accounting for the interac- 
tion between Cooper pairs occupying different quantum 
states. Interaction terms that depend on the phases of 



the coefficients appear when Li + L3 = 2L2 (|], but they 
are not considered here since these subspaces do not play 
any role in our discussion. Stationary solutions of the 
projected functional (|^) can thus be found analytically 
with respect to \Cl\ and numerically with respect to B. 
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FIG. 1. CT magnetization curves associated with local 
minima of different topological charges Q for a disk radius 
R = 55(0) and k — 1. The equilibrium magnetization is rep- 
resented by the thick solid line. Inset: Energy along the local 
minimum Q — 9 compared to the energy of the saddle point 
that separates it from the closest neighboring local minima. 

As we have anticipated, in the expression n does 
not take on the bulk nominal value for Al, but an effec- 
tive one that takes into account the geometry of the disk. 
It is very difficult to estimate this effective value, but the 
experimental evidence of the existence of multiple- vortex 
structures for the R = 5^(0) disk along with recent de- 
tailed numerical calculations ||^ indicate that k ^ 1. We 
will consider n — 1 throughout. The main results of this 
work are shown in Fig. ^ which shows the magnetization 
associated with all the energetically stable stationary so- 
lutions, i.e., all the local minima. These local minima 
are characterized by a vorticity or topological charge Q 
which is defined as the number of times that the phase 
of the order parameter winds around 27r as we complete 
one circumference moving along the inner surface of the 
disk. Different curves correspond to different values of 
Q. Along any of these curves the topological charge dis- 
tributes itself in a variety of ways. For large Q {Q > 12), 
the local minimum is always a giant vortex with L = Q. 
This local minimum is separated from neighboring lo- 
cal minima with L = Q + 1 and L = Q — 1 by saddle 
points which appear as energetically unstable stationary 
solutions in the subspaces {Q, Q -f- 1} and {Q — 1, Q}, re- 
spectively. Figure ^(b) shows an example of the modulus 
of the order parameter at the saddle point that separates 
the giant vortex L = 9 from the i = 10 close to the 
high-field end of the Q — 9 curve. Notice the strong de- 
pletion of the order parameter at an (arbitrary) point on 
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the surface [0 . Higher-energy saddle points also appear 
separating the giant vortex L — Q from other possible 
local minima L — Q + 2, L^Q — 2, etc., which fully 
guarantees its energetic stability. At the low- and high- 
field ends of the curve the local minimum L ^ Q merges 
with the saddle points {Q — 1,Q} and {Q,Q + 1}, respec- 
tively (see, for instance, the high-field side of the inset in 
Fig. |l|). In the presence of dissipation the system will 
be driven toward the neighboring local minimum. (A 
dynamical analysis |l^] of this process would be interest- 
ing, but it is beyond the scope of this work.) For smaller 
Q {Q < 12) a giant vortex can be typically found close 
to the high-field end of the curve, but, as we move along 
the curve towards lower values of H, we cross the critical 
field where the multiple vortex solution in the form of a 
ring, {0, Q}, becomes energetically favorable There 
is no barrier separating them and a weak second order 
phase transition takes place. Similarly, there are sad- 
dle points separating the solution {0, Q} from the local 
minima {0, Q—l} and {0, Q -|- 1} on neighboring curves. 
These barriers appear as stationary solutions in the sub- 
spaces {0, Q — l, Q} and {0, Q, Q-l-1}, respectively. There 
is also a barrier separating the local minimum {0, Q} 
from a local minimum {1,(5} which may become ener- 
getically favorable as we move upwards along the curve. 
There is a saddle point separating them which appears 
as a stationary solution in the subspace {0, 1, Q}. These 
saddle points are typically one or two orders of magnitude 
smaller than those separating states with different topo- 
logical charges, although, at the energy crossing, they 
still range from « 50 K (Q = 8) to « 0.1 K (Q = 11). The 
experimental temperature is 0.4 K, which means that the 
system could evolve along the state {0, Q} in the absence 
of fiuctuations. Otherwise, the system would experience 
a weak first order transition on crossing the barrier and 
changing the vortex structure into a ring with a vortex 
in the middle jl^. We have chosen this possibility in 
our plot, although these first order transitions barely re- 
flect in the magnetization. More structural changes and 
more weak first order transitions can take place as one 
moves along the curve on decreasing H. At the low- 
field end of the curve, the saddle point separating the 
local minimum {Li, . . . , Q} from {Li, . . . ,Q — 1}, i.e., 
the {Li, . . . ,Q — 1, Q} stationary solution, merges with 
the local minimum {Li, . . . ,Q} (see, for instance, the 
low-field side of the inset in Fig. |^) and the energetic 
instability sets in. Before reaching the low-field end, as 
Fig. ||(a) shows, the superconducting density exhibits 
a strong depletion at some point on the surface when 
crossing the saddle point. 

As Fig. 1 shows, the magnetization associated with the 
local minima changes sign at some intermediate value of 
H (which depends on Q). In other words, the super- 
current reverses the direction in which it flows. At this 
point the height of the vortex escape barrier approaches, 
typically, its maximum value (see inset in Fig. 1 where 



the zero-magnetization point coincides with the mini- 
mum in the free energy). As already mentioned, these 
barriers can be several orders of magnitude the experi- 
mental temperature. Thus, whatever mechanisms may 
be operating to decrease the topological charge are un- 
likely to be efficient enough to prevent the appearance of 
the paramagnetic response on decreasing H. In fact, this 
is what is observed in the experiment. The sign change 
in the response occurs approximately when the dominant 
eigenfunction in the expansion of the order parameter, 
Ln = Qi crosses the minimum of the band structure 
^l{B) and reverses its group velocity. As was shown in 
Refs. |9|,[l0|, the minimum in the band structure comes 
about due to the boundary conditions imposed to the 
components of the order parameter. Notice that detec- 
tor size effects need not be invoked Q to account for the 
paramagnetic response. 

(a) 




FIG. 2. Modulus of the order parameter at the saddle 
points separating the local minimum (a) {1,9} from {1,8} 
close the low-field end of the (5 = 9 curve and (b) {9} from 
{10} close the high-field end of the same curve. 

The experimental magnetization curves 1^ are, to a 
large extent, similar to the theoretical ones for small Q. 
However, for large where our method is expected to 
be more reliable, there are significant discrepancies. The 
derivative with respect to the field of the magnetization 
curves changes sign close to the low- and high-field ends 
of the curves. Furthermore, neighboring curves even get 
to cross each other. This is never seen in the experiment 
which seems to indicate that a vortex can escape or 
enter the disk before the saddle point or barrier disap- 
pears. Since thermal activation is only efficient very close 
to both ends of the curves, either experimental noise or 
quantum fluctuations may be ultimately responsible for 
the vorticity change. 
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Even if the energetic instability on the diamagnetic 
side is preempted by some of the relaxation processes 
mentioned above, the magnetization associated with the 
global minimum (thick line in Fig. |l]) is not likely to be 
observed for increasing field without intentional relax- 
ation. It has been suggested in the literature that sur- 
face roughness is responsible for the destruction of the 
saddle points associated to the surface barrier. The sad- 
dle points preventing the escape or entrance of vortices 
have the same origin and surface roughness should affect 
them similarly. Thus, in our view, there are no reasons 
for the system to follow the ground state on decreasing 
H and it is expected to continue along a given curve un- 
til the onset of the energetic instability or until vortex 
entrance rates increase to typical measurement times. 
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FIG. 3. (a) FC magnetization for a disk of radius R — 5^(0) 
and K = 1 at different values of the external magnetic field. 
Each curve corresponds to the topological charge Qn of the 
giant vortex that nucleates at the highest critical temperature, 
(b) The same as in (a), but for a topological charge Q„ + 1. 

Finally, we would like to comment on the FC mea- 
surements. The current understanding of the FC re- 
sults is summarized in a work by Moshchalkov et al. fl^ ] 
which attributes the FC paramagnetic response to a flux- 
compression phenomenon. Figure ^ shows the magneti- 
zation as a function of temperature for different values of 
the magnetic field. The usual phcnomcnological temper- 
ature scaling of the parameters in the Ginzburg-Landau 
functional (|]) has been considered [|^. Each curve in 
Fig. ^(a) corresponds to the topological charge Q„ of the 
giant vortex that nucleates at the highest possible criti- 
cal temperature for each chosen field. Qn is maintained 
along the different curves down to T = due to presence 
of the energy barriers discussed above which prevent the 
change of vorticity. The response is always diamagnetic 
in clear contrast with the FC data. Figure |^(b) shows 
magnetization curves for a topological charge Qn + 1. 
Alternating paramagnetic and diamagnetic behaviors are 
obtained as a function of H and a low-temperature satu- 
ration of the paramagnetic response occurs due to explo- 
sion of the giant vortex as suggested by Moshchalkov et 
al. pj]. This behavior is in remarkable agreement with 



the FC data ||^ which seems to suggest that either ther- 
mal fluctuations close to the critical temperature or sur- 
face roughness favor the nucleation of giant vortices with 
a higher topological charge than that expected from plain 
Ginzburg-Landau theory [ p^ . Notice again that, in our 
approach, the magnetic field is uniform in space which 
suggests that flux compression is not essential as 

far as the existence of paramagnetism is concerned. Still, 
the origin of the PME in the FC measurements remains 
an open issue. 

In closing, we have addressed some fundamental ques- 
tions posed by the experiment of Geim et al. |3|. Sad- 
dle points or barriers of the Ginzburg-Landau functional 
have been found based on a projection technique. This 
has allowed us to obtain metastable magnetization curves 
and to gain insight into the controversial paramagnetic 
response both in FC and CT magnetization measure- 
ments. 
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